ON ASYMPTOTIC STABILITY OF NONCHARACTERISTIC VISCOUS 

BOUNDARY LAYERS 



TOAN NGUYEN 



Abstract. We extend our recent work with K. Zumbrun on long-time stability of multi- 
dimensional noncharacteristic viscous boundary layers of a class of symmetrizable hyperbolic- 
parabolic systems. Our main improvements are (i) to establish the stability for a larger 
class of systems in dimensions d > 2, yielding the result for certain magnetohydrodynamics 
(MHD) layers; (ii) to drop a technical assumption on the so-called glancing set which was 
required in previous works. We also provide a different proof of low-frequency estimates 
by employing the method of Kreiss' symmetrizers, replacing the one relying on detailed 
derivation of pointwise bounds on the resolvent kernel. 
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1. Introduction 



Boundary layers occur in many physical settings, such as gas dynamics and magneto- 
hydrodynamics (MHD) with inflow or outflow boundary conditions, for example the flow 
around an airfoil with micro-suction or blowing. Layers satisfying such boundary conditions 
are called noncharacteristic layers; see, for example, the physical discussion in [S| ISGKQ] . 
See also jGMWZSl lYZl INZTI lNZ2l [Z5] for further discussion. 

In this paper, we study the stability of boundary layers assuming that the layer is non- 
characteristic. Specifically, we consider a boundary layer, or stationary solution, connecting 
the endstate ?7+: 

(1.1) U = U{xi), lim U{xi) = U+. 

of a general system of viscous conservation laws on the quarter-space 

(1.2) Ut + ^F^{U)x^=Y,iB'''{U)U.,)x,, x€Mt t > 0, 

j jk 

U,F^ G M", B^^ G M"-^", with initial data jj{x,0) = Uo{x) and boundary conditions as 
specified in (B) below. 

An fundamental question is to establish asymptotic stability of these solutions under 
perturbation of the initial or boundary data. This question has been investigated in \GR\ 
IMeZli IGMWZSi IGMWZ61 lYZl iNZTl |NZ2 ^ for arbitrary-amplitude boundary-layers using 
Evans function techniques, with the result that linearized and nonlinear stability reduce to 
a generalized spectral stability, or Evans stability, condition. See also the small-amplitude 
results of [GGl lR3l [MM iKNZl iKaK] obtained by energy methods. 

In the current paper, as in [Nl| for the shock cases, we apply the method of Kreiss' 
symmetrizers to provide a different proof of estimates on low-frequency part of the solution 
operator, which allows us to extend the existing stability result in [NZ2j to a larger class 
of symmetrizable systems including MHD equations, yielding the result for certain MHD 
layers. We are also able to drop a technical assumption (H4) that was required in previous 
analysis of [ZH [Z3l [M IGMWZll [NZ2] . 



1.1. Equations and assumptions. We consider the general hyperbolic-parabolic system 
of conservation laws ()1.2p in conserved variable U, with 



ull I ^ ^ - \ dk iJk 



6f 6iV' 

~I ^ ^n-r^ ^11 ^ ^r^ 

jk 

Following |MaZ3l [MaZ4l IZ31 IZ4j . we assume that equations p.2p can be written, alter- 
natively, after a triangular change of coordinates 

(1.3) W := W(0) = U|^'ln^ , 
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in the quasilinear, partially symmetric hyperbolic-parabolic form 
(1.4) A^Wt + A^W., = Y.iB^'W.,)., + G, 

3 jk 

where, defining W+ := W{U+), 

(Al) {W+) , , are symmetric, A^ block diagonal, > 6*0 > 0, 

(A2) for each ^ G M'^ \ {0}, no eigenvector of CjA^ iA°y^(W+) lies in the kernel of 

(A3) W'' = E^%-^;t > and G = (^^ with g{W,,W,) = 0{\W,\^). 

Along with the above structural assumptions, we make the following technical hypotheses: 

(HO) F^,B^'',A°,A^,&'',W{-),g{-,-) G C^ with s > [{d - l)/2] + 5 in our analysis of 
linearized stability, and s > s{d) := [{d — l)/2] + 7 in our analysis of nonlinear stability. 

(HI) A\^ is either strictly positive or strictly negative, that is, either A\^ > > 0, or 
A\^ < —6*1 < 0. (We shall call these cases the inflow case or outflow case, correspondingly.) 

(H2) The eigenvalues of dF^{U+) are distinct and nonzero. 

(H3) The eigenvalues of Y^j ^jdF^ {Uj^) have constant multiplicity with respect to ^ G M'^, 

Alternative Hypothesis H3'. The constant multiplicity condition in Hypothesis (H3) 
holds for the compressible Navier Stokes equations whenever is hyperbolic. We are able 
to treat symmetric dissipative systems like the equations of viscous MHD, for which the 
constant multiplicity condition fails, under the following relaxed hypothesis. 

(H3') The eigenvalues of S^jdF^ (U^) are either semisimple and of constant multiplicity 
or totally nonglancing in the sense of [GMWZ6] . Definition 4.3. 

Additional Hypothesis H4' (in 3D). In the treatment of the three-dimensional case, 
the analysis turns out to be quite delicate and we are able to establish the stability under the 
following additional (generic) hypothesis (see Remark 13.41 and Appendix |A] for discussions 
of this condition): 

(H4') In the case the eigenvalue Xk{C) of J2j ^j'^F'^ {U+) is semisimple and of constant 
multiplicity, we assume further that V^A^ ^ when d^-^^Xk = 0, ^ 7^ 0. 

Remark 1.1. Here we stress that we are able to drop the following structural assumption, 
which is needed for the earlier analyses of \Z2\ \Z3\ IZ41 INZ2] . 

(H4) The set of branch points of the eigenvalues of {A^)~'^{iTA^ + r G M, 

^ G M'^"^ is the (possibly intersecting) union of finitely many smooth curves r = on 
which the branching eigenvalue has constant multiplicity Sq (by definition > 2). 
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We also assume: 

(B) Dirichlet boundary conditions in VF-coordinates: 

(1.5) {w^,w^^){0, X, t) = h{x, t) := {hi, h2){x, t) 
for the inflow case, and 

(1.6) w^\Q,x,t) = h{x,t) 
for the outflow case, with x = {xi,x) G W^. 

This is sufficient for the main physical applications; the situation of more general, Neu- 
mann and mixed-type boundary conditions on the parabolic variable w^^ can be treated as 
discussed in |GMWZ5l[G"MWZ6j . 

1.2. The Evans condition and strong spectral stability. A necessary condition for 
linearized stability is weak spectral stability, defined as nonexistence of unstable spectra 
KA > of the linearized operator L about the wave. As described in \L2\ IZ3] . this is 
equivalent to nonvanishing for all ^ G M'^"^, 5RA > of the Evans function 

DlH A), 

a Wronskian associated with the family of eigenvalue ODE obtained by Fourier transform in 
the directions x:={x2,..., xj). See [Z3l [GMWZ51 IC^MWZGl lNZ2] for further discussion. 

Definition 1.2. We define strong spectral stability as uniform Evans stability: 

(D) \DL{l\)\>e{C)>{) 

for (I, A) on bounded subsets C C {| G W^-^ , KA > 0} \ {0}. 

For the class of equations we consider, this is equivalent to the uniform Evans condition of 
[GMWZSl IGMWZ6] , which includes an additional high-frequency condition that for these 
equations is always satisfied (see Proposition 3.8, |GMWZ5] ). A fundamental result proved 
in [GMWZ5] is that small-amplitude noncharacteristic boundary-layers are always strongly 
spectrally stable. 

Proposition 1.3 ([GMWZ5]). Assuming (A1)-(A3), (H0)-(H2), (H3'), (B) for some fixed 
endstate (or compact set of endstates) U^, boundary layers with amplitude 

\\U — ^^+||l°o[0,+oo] 
sufficiently small satisfy the strong spectral stability condition (D). 



As demonstrated in [SZl IZ5j . stability of large-amplitude boundary layers may fail for 
the class of equations considered here, even in a single space dimension, so there is no 
such general theorem in the large-amplitude case. Stability of large-amplitude boundary- 
layers may be checked efficiently by numerical Evans computations; see, e.g., [HLZl [CHNZl 
HLyZltlHL^ . 
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1.3. Main results. Our main results are as follows. 

Theorem 1.4 (Linearized stability). Assuming (A1)-(A3), (H0)-(H2), (H3'), (H4'), (B), 
and (D), we obtain the asymptotic H stability in dimensions d > 3, 

and any 2 < p < oo, with rates of decay 

^^^^ ![/(i)!L2 <c(i + t)-^-1t/o|Lw, 

for some e > 0, provided that the initial perturbations Uq are in Li n //[(<i-i)/2]+2^ and zero 
boundary perturbations. 

Theorem 1.5 (Nonlinear stability). Assuming (Al)-(A3), (H0)-(H2), (H3'), (H/), (B), 
and (D), we obtain the asymptotic L} R — > n stability in dimensions d > 3, for 
s > s{d) as defined in (HO), and any 2 < p < oo, with rates of decay 

\U{t) -U\l.< C{1 + tr'^^'~'/^^''^.-'\Uo\L^nHs 

\U{t) - U\hs < C{1 + t)-^-'\Uo\LinHs, 

for some e > 0, provided that the initial perturbations Uq := Uq — U are sufficiently small 
in n and zero boundary perturbations. 

Remark 1.6. As will be seen in the proof, the assumption (H4') can be dropped in the 
case (i > 4, though we then lose the factor in the decay rate. 

Our final main result gives the stability for the two-dimensional case that is not covered 
by the above theorems. We remark here that as shown in [Z2l IZ3] . Hypothesis (H4) is 
automatically satisfied in dimensions d = 1 , 2 and in any dimension for rotationally invariant 
problems. Thus, in treating the two-dimensional case, we assume this hypothesis without 
making any further restriction on structure of the systems. Also since the proof does not 
depend on dimension d, we state the theorem in a general form as follows. 

Theorem 1.7 (Two-dimensional case or cases with (H4)). Assume (A1)-(A3), (H0)-(H2), 
(H3f), (H4), (B), and (D). We obtain asymptotic n H"" ^ LP f] H"" stability of U as a 
solution of (|1.2p in dimension d >2, for s > s{d) as defined in (HO), and any 2 < p < oo, 
with rates of decay 

\Uit) - U\l. < C{1 + tyl^'-'M+'/'PlUQlLinH" 

\U{t)-U\Hs <Cil+t)-^\UQ\LinHs, 

provided that the initial perturbations Uq := Uq — U are sufficiently small in D H^ and 
zero boundary perturbations. Similar statement holds for linearized stability. 

Remark 1.8. The same results can be also obtained for nonzero boundary perturbations 
as treated in |NZ2| . In fact, in |NZ2] . though a bit of tricky, it has been already shown that 
estimates on solution operator (see Proposition 12. ip for homogenous boundary conditions 
are enough to treat nonzero boundary perturbations. Thus for sake of simplicity, we only 
treat zero boundary perturbations in the current paper. 
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Combining Theorems 11.4^ 11.51 11-71 and Proposition 11.31 we obtain the following small- 
amplitude stability result. 

Corollary 1.9. Assuming (Al)-(A3), (H0)-(H2), (E3), (B) for some fixed endstate (or 
compact set of endstate s) boundary layers with amplitude 

sufficiently small are linearly and nonlinearly stable in the sense of Theorems I j.5|. and 

[13 

2. Nonlinear stability 
The linearized equations of ()1.2p about the profile U are 



(2.1) Ut = LU := Y.iB^^U^J,^ - ^(A^U),^ 

j,k j 

with initial data ?7(0) = Uq. Then, we obtain the following proposition, extending Propo- 
sition 3.5 of [NZ2] under our weaker assumptions. 

Proposition 2.1. Under the hypotheses of Theorem \1.5l the solution operator S{t) := e^* 
of the linearized equations may be decomposed into low frequency and high frequency parts 
(see below) as S{t) = Si{t) + S2{t) satisfying 

\Sl{t)d^.ld§f\,. <C{1 + t)-(^-2)/4-e/2-|,|/2|^|^^ ^ c{l + t)-^''-^y'-^l^\fy^^^ 

(2.2) |5i(t)af^af /i^.,. <c{i + t)-('^-i)/4-/2H/3|/2|^| ^ ^ 

\s,m\dlf\L- <c{i + t)~(<^-i)/w2-i/3i/2|^|^^ + c{i + t)-^'^-^y^-^/^\f\,.^ 



x.x-^ 



for some e > and 13 = /3) with f3i = 0, 1, and 

(2.3) \dl\dlS2{t)f\L2 < Ce-^^*|/l^/hil+m+3, 

for 7 = (71, 7) with 71 = 0, 1. 

We shall give a proof of Proposition 12.11 in Section [3l For the rest of this section, we give 
a rather straightforward proof of the first two main theorems using estimates of the solution 
operator stated in Proposition 12. H following nonlinear arguments of |Z3t INZ2] . 

2.1. Proof of linearized stability. Applying estimates on low- and high-frequency oper- 
ators Si{t) and S2{t) obtained in Proposition 12.11 we obtain 

\U{t)\L2 < \Sl{t)Uo\L2 + \S2{t)UQ\L2 

(2.4) < C{l+t)-'-^-i[\Uo\L^ + |t/oLi,oo ] + Ce-'^'\Uo\m 

< C{l+t)-^-^UQ\Llr^H^ 
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and (together with Sobolev embedding) 

\U{t)\L^ < \Si{t)Uo\L^ + \S2{t)Uo\L^ 

< C7(l + t)-^'f [|[/o[ii + |C/oLi,oc] +C|52(t)C/olHM-i)/2)+2 

(2.5) , 

< C{l + t)-—'2[\Uo\Li + |C/olii,-] + Ce-^''*|?7o|^,[(.-i)/2,+2 

d~l e 

<C{l+t) 2 2\Uo\^i^j^[(d-l)/2]+2. 

These prove the bounds as stated in the theorem for p = 2 and p = oo. For 2 < p < oo, we 
use the interpolation inequahty between and 

2.2. Proof of nonlinear stability. Defining the perturbation variable U := U — U , we 
obtain the nonlinear perturbation equations 

(2.6) Ut-LU = Y,QHU,U,).,, 

3 

where 

Q^{U,U^) = 0{\U\\U^\ + \U\'') 

(2.7) Q'{U,U,\^=0{\U\\U^\ + \U\\U^^\ + \U^\'') 

Q'iU,U,),^,, = 0{\U\\U,,\ + \U^\\U,,\ + + \U\\U,,,\) 

so long as remains bounded. 

Applying the Duhamel principle to (j2.6p . we obtain 

(2.8) U{x,t) =S{t)Uo + S{t - s)Y,d.,QHU,U^)ds 

where U{x,0) = Uo{x). 
Proof of Theorem \1.5[ Define 

at) := sup(|^7(.)|^2(l + 3)"^+^ + \U{s)\ls^{1 + s)^+^ 

(2.9) 

+ (|C/(5)| + |[/,(S)|)^2,^ (1 + .)^ 



X ,x 



We shall prove here that for all t > for which a solution exists with C(t) uniformly 
bounded by some fixed, sufficiently small constant, there holds 

(2.10) C{t)<C{\Uo\L^nHs+m^)- 
This bound together with continuity of C{t) implies that 

(2.11) at) < 2C\Uo\L^nHs 

for t > 0, provided that iC/olLin/f^ < 1/4C^. This would complete the proof of the bounds 
as claimed in the theorem, and thus give the main theorem. 

By standard short-time theory/local well-posedness in H^, and the standard principle 
of continuation, there exists a solution U S on the open time- interval for which \U\h<' 
remains bounded, and on this interval C(i) is well-defined and continuous. Now, let [0,T) 
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be the maximal interval on which |{7|/fs remains strictly bounded by some fixed, sufficiently 
small constant 5 > 0. Recalling the following energy estimate (see Proposition 4.1 of |NZ2j ) 
and the Sobolev embeding inequality |[/|^^2,oo < C|C/|/fs, we have 

\U{t)\h < Ce-'^\Uo\h+C f e-'^'-^\U{T)\l,dT 

(2.12) Jo 

< c{\u^\is + at?){i + 1)-('^-2)/2-2.. 

and so the solution continues so long as C, remains small, with bound (12. lip , yielding exis- 
tence and the claimed bounds. 

Thus, it remains to prove the claim (I2.10p . First by ()2.8p . we obtain 

\U{t)\L2 <\S{t)Uo\L^+ f \Si{t-s)d^^Q^{s)\L2ds 

(2.13) ^ •'^ 

+ / \S2{t-s)d^^Q^{s)\L2ds 
Jo 

where \S{t)Uo\L2 < C(l + t)"^~'|C/o|iin//3 and 
\Si{t-s)d,^Q^{s)\L2ds 

<C {l + t- s)-^-"2--\Q^{s)\l. + (1 + s)-^-'\Q\s)l.,o. ds 

Jo ""'^i 
<C {l + t-s)-'^-^2-'\U\'jJ^ + {l + t-s)-'^-'(\U\\,^ +\U^\^2,o.)ds 

Jo ^ i.a^i 4,11 / 

< C{\Uo\h + atf) r \{l + t- s)-^-5-(l + s)-^-2^ 



+ (l+t-s)-— "^(1 + s)- 



d-1 



2 2. 



ds 



d~2 



<c{i+t)--^-^{\uo\is+at?' 



and 



\S2{t-s)d^^Q^{s)\L2ds 

< f e-'^'-'\d,^Q\s)\H^ds 
Jo 

<C [ e-^^'-''^\U\lsds 
Jo 

< C{\Uo\h + C{tf) f e-^(*-^)(l + sy'-^-^^ds 

Jo 



<c{i+t)-^-^%\uo\i.+at?)- 

Therefore, combining these above estimates yields 
(2.14) |C/(t)U2(l + t)^+^ < C{\Uo\L^nHs+Q{tf ). 
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Similarly, we can obtain estimates for other norms of U in definition of C, and finish the 
proof of claim ()2.10p and thus the main theorem. □ 

Remark 2.2. The decaying factor is crucial in above analysis when d = 3. In fact, the 
main difficulty here comparing with the shock cases in |N1] is to obtain a refined bound of 
solutions in L°°. See further discussion in Section [3] below. 

3. Linearized estimates 

In this section, we shall give a proof of Proposition 12.11 or bounds on Si{t) and 52 (t), 
where we use the same decomposition of solution operator S{t) = Si{t)+S2{t) as in 



3.1. High— frequency estimate. We first observe that our relaxed Hypothesis (113') and 
the dropped Hypothesis (H4) only play a role in low-frequency regimes. Thus, in course of 
obtaining the high-frequency estimate (j2.3p . we make here the same assumptions as were 
made in [NZ2] , and therefore the same estimate remains valid as claimed in (j2.3p under our 
current assumptions. We omit to repeat its proof here, and refer the reader to the paper 
[NZ2] . Proposition 3.6. 

In the remaining of this section, we shall focus on proving the bounds on low-frequency 
part Si{t) of linearized solution operator. 

Taking the Fourier transform in a; := (x2, . . . , x^) of linearized equation (j2.ip . we obtain 
a family of eigenvalue ODE 

LoU 

XU = Lfl := {BiiU'y - {AiUy -i^Aj^jU + iJ2Bji^jU' 

(3.1) m m 

3.2. The GMWZ's stability estimate. Let U = {u^ ,u^^)'^ be a solution of resolvent 
equation (L^ — X)U = f. Following [Z3, GMWZ6| . consider the variable W as usual 

/ V'^ 

W : = 

with := A^u^,w^^ := b\^u^ + bl^u^^ , A^ := A\^ - ^l2(^2^)~^^P- Then we can write 
equations of as a first order system 

.3 2^ d,,W = gixi,X,i)W + F 

rP^ = on xi = 0. 

For small or bounded frequencies (A,^), we use the MZ conjugation lemma (see [MeZl^ 
IMeZ3] ). That is, given any (A, ^) G 1^"^+^, there is a smooth invertible matrix <I?(xi, A, ^) for 
xi > and (A,.^) in a small neighborhood of (A, .^), such that (|3.2p is equivalent to 

(3.3) d.,,Y = g+{X,i)Y + F, t{X,i)Y = 

where g+{X, 1) := §{+00, A, 1), W = <^Y,F = <^'^F and FY := T^F. 
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Next, there are smooth matrices ^(A, such that 

(3.4) v-'g+v ° 





with blocks H{X,^) and P(A,^) satisfying the eigenvalues /x of P in {|3?e/i| > c > 0} and 

d 

Ho{X,0 ■■ = -{Al)-'{{iT + j)Al + Y,i^,A^), 

i=2 

with A = 7 + ir. 

Define variables Z = {uH,upf as W = <^Y = ^VZ, fZ := T^VZ, and {/hJpV = 
V-^F. We have 

(3-5) 5.,M = f^ DM + rz = o. 



Up J I Pj \up J \fp 

Then the maximal stability estimate for the low frequency regimes in [GMWZ6] states 
that 

(3.6) (7 + p^)\uH\l2 + \up\l, + lunm'' + l«p(0)P < (Ifnl \uh\) + (|/p|, \up\). 

We note that in the final step there in [GMWZT] . the standard Young's inequality has 
been used to absorb all terms of {uh,up) into the left-hand side, leaving the norm of F 
alone in the right hand side. For our purpose, we shall keep it as stated in (j3.6p . Here, by 
/ < 5(, we mean f < Cg, foi some C independent of parameter p. 

We remark also that as shown in [GMWZT] , all of coordinate transformation matrices are 
uniformly bounded. Thus a bound on Z = {uh,up)'^ would yield a corresponding bound 
on the solution U. 

3.3. and L°° resolvent bounds. Changing variables as above and taking the inner 
product of each equation in ()3.5p against uh and up, respectively, and integrating the 
results over [0,xi], for xi > 0, we obtain 

1 1 f^^ 

-\uh{xi)\'^ = -\uh{0)\'^ + '!Re / {H{X,^)UH ■ UH + Ih ■ UH)dz, 

(3-7) ^ ^ 

-\up{xi)\^ = -\up{0)\^ + / {P{X,i)up-up + fp-up)dz. 
^ I Jo 

This together with the facts that \H\ < C p and \P\ <C yields 

\uh\\c^(^-,)<\uh{^)?' + P\uh\\2 + {\fH\AuH\), 

(3.8) 2 2 2 

kp|L-{a'i) ^ hp(0)l + \up\l2 + (|/p|, |np|), 
and thus in view of (j3.6p gives 

(3.9) (7 + p^)\uh\\2 + \up\\2 + p\uh\\oc. + \up\\oo < {\fH\, \uh\) + {\fp\, \up\). 
Now applying the Young's inequality, we get 

ilfnl, \uh\) + (I/p|, \up\) < {e\up\lc^ + Ce\fp\li) + [ep\uH\l^ + "^l/^^lii 
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and thus for e sufficiently small, together with (j3.9 

(3.10) ' "2m..__|2 , L._|2 , „L,__|2 

Therefore in term of Z = {uh,upY, 

(3.11) |Z|icx,(^j) < Cp'VIli and 



(7 + P'')\uH\'i2 + \up\l2 + p\uh\Io. + \up\l^ < ^Ifnlli + \fp\li- 



Unfortunately, unlike the shock cases (see [Nlj ). bounds (|3.1ip are not enough for our 
need to close the analysis in dimension d = 3. See Remark 12.21 In the following subsection, 
we shall derive better bounds for Z in both L°° and norms. 

3.4. Refined and L°° resolvent bounds. With the same notations as above, we prove 
in this subsection that there hold refined resolvent bounds: 

(3.12) < p-i+^d/ki + l/U^) and \Z\l2^,^) < p-'/^+'{\f\L^ + |/|l-) 

for some small e > 0. We stress here that a refined factor p'' in is crucial in our analysis 
for three-dimensional case. See Remark 12.21 

Assumption (H3') implies the following block structure (see |MeZ3l [GMWZ6] ) . Here, we 

use the polar coordinate notation ( = (r, 7,^),^ = pC where ( = (f,7,^) and ( € S'^. 

Proposition 3.1 (Block structure; |GMWZ6] ). For all C with 7 > there is a neighborhood 

(jj of (C, 0) in 5"^ X M+ and there are C°° matrices T{(, p) on u) such that T~^HqT has the 
block diagonal structure 

(3.13) T-^HoT = Hb{C, p) = pHsiC, p) 
with 

" Qi 

(3.14) Hb{Lp)= ••• (C,P) 

Qp 

with diagonal blocks of size ^ such that: 

(i) (Elliptic modes) ^Qk is either positive definite or negative definite. 

(a) (Hyperbolic modes) = 1, ^Qk = when 7 = /) = 0, and dx^{^Qk)dp(JSiQk) > 0. 

(Hi) (Glancing modes) Vk > 1? Qk has the following form: 



(3.15) 

where a :- 

(3.16) 



Qfc(C, p) = z(/i Jd + J) + iaQ'fcd) + 0(7 + p). 



11 -f I, 



J 



1 





Q'kCi) ■■- 



(12 













7^ 0, and the lower left hand corner a of Qk satisfies dA^{^a)dp{^a) > 0. 
(iv) (Totally nonglancing modes) > 1, eigenvalue of Qk, when j = p = 0, is totally 
nonglancing, see Definition 4-3, jGMWZB] . 
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Proof. For a proof, see for example [Met], Theorem 8.3.1. It is also straightforward to see 
that for the case (iii), 

Q,,(|) = lV^-Z)fc(C,i,)|=c|V^-Afc(0|, 

where c is a nonzero constant, -Dfc(C)Ci) is defined as det{iQk{C) +^ild), and XkiC) is the 
zero of -Da;(C)^i) (recalling ( = (A,^)) satisfying 

%Afc = ... = a^^-^Afc = 0, d'^'^Xk + at ilQ. 

Thus, assumption (H4') guarantees the nonvanishing of q^^,. We skip the proof of other 
facts. □ 

We shall treat each mode in turn. The following simple lemma may be found useful. 

Lemma 3.2. Let U he a solution of dzU = QU + F with U{+oo) = 0. Assume that there 
is a positive [resp., negative] symmetric matrix S such that 

(3.17) ^SQ :=^{SQ + Q*S*)>0Id 

for some 9 > 0, and S > Id [resp., —S > Id]. Then there holds 

, , ^ <\F\\, 

[resp., \U\l^ + 9\U\l, < [[/(O)^ + \F\l^ ]. 

Proof. Taking the inner product of the equation of U against SU and integrating the result 
over [xi, oo] for the first case [resp., [0, xi] for the second case], we easily obtain the lemma. 

□ 

Thanks to Proposition 13. 11 we can decompose U as follows 

(3.19) U = up + UH, + + UHg + UHt , 

corresponding to parabolic, elliptic, hyperbolic, glancing, or totally nonglancing modes. 

3.4.1. Parabolic modes. Since spectrum of P is away from the imaginary axis, we can assume 
that 



P+ 

P_ 



with ib9?P± > c > 0. Therefore applying Lemma 13.21 with S = Id or —Id yields 

I |2 I I |2 ^ I 17 |2 

|up_|ioc + \up_\l2 < \up_{0)\^ + \Fp_\li. 
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3.4.2. Elliptic modes. This is case (i) in Proposition 13.11 when the spectrum of Qk hes in 

{3?e/i > 6} [resp., {^efi < -6}]. 

In this case, there are positive symmetric matrices S''^((",/j), C°° on a neighborhood uj of 
(C, 0) and such that 

^S''Q^ > cid [resp., - ^S^Q^ > cid] 
for c > 0. Thus, Lemma 13.21 again yields 

\UH,_\1^ +P\UH,_\12 < \uH,_iO)\'^ + \Fh,_\Ii. 

3.4.3. Hyperbolic modes. This is case (ii) in Proposition [3?T1 In this shown in |Metj 
we can write 

(3.22) Q\Cp)=qHC)Id + pn''{C,p) 

where is purely imaginary when 7 = 0, g'^ := d^^eq^{(^) does not vanish, and the 
spectrum of q^lZ^iCjO) is contained in the half space {^ep > 0}. Therefore, when q'' > 
[resp., q^ < 0] and thus for (C>7) sufficiently close to (C,0) 

^eq'^ > C7, [resp., ^eq^ < -cj], 
we have positive symmetric matrices S'^{(,p) satisfying 

^S'^Q'' > c(7 + p)Id [resp., - ^S^Q^ > c{j + p)Id] 
for c > 0. Thus, again by Lemma 13.21 we obtain 
^^23) K^^h+li- + {l + P^)\^Hn+\l2 < \Fh^+\Iu 

Wh^M^ + in + p')\uH^-\h ^ Wh,A^)? + \Fh^M.. 

3.4.4. Totally nonglancing modes. This is case (iv) in Proposition 13.11 As constructed in 
[GMWZH] ■ there exist symmetrizers S'^ that are positive [resp. negative] definite when the 
mode is totally incoming [resp. outgoing]. Denote UHt+ [resp., UHt-] associated with totally 
incoming [resp. outgoing] modes. Then similarly as in above, we also have 



'uHtM^ + {l + P^)\uHtM2 < \uH,^m^ + \Fm.\ll- 



Thus, putting these estimates together with noting that the stability estimate (j3.6p al- 
ready gives a bound on |u(0)|, we easily obtain sharp bounds on u in and for all 
above cases: 

(3-25) Wkll^ + P^Wk\\2 < \ f\\i + \uHg\L'^\f\L^, 

for allk = P,He,Hh,Ht. 
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3.4.5. Glancing modes. Hence, we remain to consider the final case: case (iii) in Proposition 
EH Recall ([3:15]) 

(3.26) QkiC, p) = i(^,Id + J) + iaQUi) + ^(7 + p) 
on a neighborhood of (C, 0), where a = |^ — ^|. We consider two cases. 

Case a. o p^ for some small e > 0. Recall that we consider the reduced system: 

(3.27) dn^^Uk = pQkiC, P)uk + fk 

with QkiC^p) having a form as in (j3.26p . It is clear that the norm of Uk remains 
unchanged under the transformation Uk to Ufce~*-fc^\ Thus, we can assume that = 0. 
Note that we have the following bounds by (j3.1ip 

(3-28) \uk\L°-(xi) < p'-^lfli^ and \uk\L\xi) ^ P'^^^IIIl^- 

To prove the refined bounds (j3.12p . we first observe that 

\dxiUk\L^ < p\uk\L°° + \fk\L^ ^ + I/Il'^j 

where the last inequality is due to (j3.28p . Now, write Uk = (wfc,i, ■ ■ ■ , Uk,uk)- Thanks to the 
special form of Qk in ()3.26p . we have 

(3.29) dx^Uk,Uk = ip<^Q'k{l)uk + Oin + p^)uk + fk- 

Taking inner product of the equation (13.29P against dxiUk^u^^ easily obtain by applying 
the standard Young's inequality: 

(3.30) \dx,Uk,uJh ^ P'^^>k\l2 + \fk\LAdx,Uk,.u\L^ < p-'^-^Vlh + l/li-- 
Similarly, for u^^y^-i satisfying 

dxiUk,Uk-i = ip(^Q'kii)uk + ^puk,,y^: + 0(7 + p'^)uk + fk, 

we have 

(3.31) \dxiUk,Uk~i\L2 < p'^^'^^Ukl'ii + p\ < Uk,Uk,dxiUk,Uk-i > \ + IfklL'^ldxiUk^UklL^^- 
Here, integration by parts and Young's inequality yield 

P\ < Uk,Uk,dxiUk,uk-l > I ^ P\dxiUk,uk\L'2\Uk,Uk-l\L2 +P|'Wfc(0)|^. 

Thus, using the refined bound (j3.30p and noting that 

l^^fc(0)P <\<f,Uk>\< l/lLikfcU- < /o^Vlii, 

we obtain 

p\ < uk,u,,dx,uk,u,-i > I < + l/Ii-) 

Therefore, applying this estimate into (j3.3ip . we get 

(3.32) \dx,nk,u,-i\l2 < p~'^mh + l/li-)- 

Using this refined bound, we can estimate the same for Uk,u^~2-, Uk,u^~3, and so on. Thus, 
we obtain a refined bound for Uk- 

(3.33) \dx,Uk\h<p-'+H\f\l + \f\l^) 
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where e may be changed in each step and smaller than the original one. This and the 
standard Sobolev imbedding yield 



(3.34) 



which proves the L°° refined bound in (13.12P for Z. Using (j3.34p into (13. 9|) . we also obtain 
the refined bound in L? as claimed in p.l2p : 



(3.35) 

for some e > 0. 



<p-^^\\f\h + \f\l^), 



Case b. o ^ p'^ for some small e in (0,1/2). We shall diagonalize this block. Recall 
that 



(3.36) 

Following 



Qk{C,p) = ^/ffcld + i 



10 

■•. 

■•. ■•. 1 





+ 0{a) 



L (Jq^k 

IGMWZT] . we diagonalize this glancing block by 



where UHg '■= UHg+ + uh^_. Here uh^^ are defined as the projections of onto the 
growing (resp. decaying) eigenspaces of QkiC^p) (I3.36P . We recall the following whose 
proof can be found in \Z2\ IZ3j or Lemma 12.1, [ GMWZT] . 



Lemma 3.3 (Lemma 12.1, [GMWZTj ). The diagonalizing transformation Tug may he cho- 
sen so that 



(3.37) \THg\<C, \ThI\<C^, 

where a, /? are defined as 



(3.38) 

and T 



fi := a 



a := a 



\T-\ I < Ca 



(l-[{z.fc+l)/2])/i.fe 



-1 



denotes the restriction of Tj^^ to subspace Hg-. In particular, /3a ^ > 1 



Simple calculations show that eigenvalues of Qfc are 

(3.39) Okj = ipk + T^k,j + o{a^^'''^), j = 0, 1, s - 1. 

Here, vr^j- = eH{q^^,ay^'^'' , with e = \ . We can further change of coordinates if necessary 
to assume that 



(3.40) 
with 

(3.41) 



Q'k ■= '^HgQkTHg = diag(afc,i, ■ • • , ak,i, ak,i+i, ■■■ , ak,Uk) 

-3?e ak,j > 0, j = 1,...,/, 

Okj > 0, j = 1 + l,-.-,i^k- 
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Hence, applying Lemma 13.21 to equations of u'^^ with S = Id or S = —Id, we easily 
obtain 

luHg+ll^ +/9min|5fte afcj||'u'^^+li2 < IF'h^+IIi, 

(3.42) 

The diagonalized boundary condition V := TaTug- By computing, we observe that 

c'-Vff I , , 

\T'u'hJ = \ruH,_\ > C-'\UH,_\ > 1 > C-'a-'\\u'Hj. 

Thus, 

(3.43) < Ca\r'u'H^_ \ < Ca(|rV| + |rV+|) < Ca\u'^\. 
Using this estimate, (|3.37p . and (|3.25p . the estimate (|3.42p yields 

(3.44) a'^hnjioo +pa"^mm|Ke akj\\uHjl2 </3^|/lii- 

Recalling that a, (3 are defined as in (j3.38p and the fact that we are in the case oi a > p'' 
for some small e > 0, we get 

(3.45) IuhJl- < C7a/3|/Ui < Cp-^^l/Ui, 
from which we obtain the refined bounds (j3.12p for this case as well. 

Remark 3.4. In case b) above, we use the nonvanishing of q^f. to make sure that (jqu^, is 
order of a in the neighborhood u; of ((", 0) so that the lower left hand entry of Qk dominates 
and thus we can be sure to diagonalize the block. Otherwise, the other entries of Qk in 
()3.36p may dominate and the behavior is not clear. The nonvanishing of q^f, is guaranteed 
by our additional Hypothesis (H4') as shown in the proof of Proposition 13.11 This is only 
place in the paper where the assumption (H4') is used. 

3.5. —> estimates. We establish the IP resolvent bounds for low frequency 

regime, restricting our attention to the surface 

(3.46) r«:={A : 3^^eA = -6'i(|||2 + |5^mAp)}, 

for 9i > Q. Taking 6i to be sufficiently small such that all earlier resolvent estimates are 
still valid on F^, with p := |(^, A)| being sufficiently small. Thus, we obtain the following: 

Proposition 3.5 (Low-frequency bounds). Under the hypotheses of Theorem \1.5l for A € 
and p := |(^, A)|, di sufficiently small, there holds the resolvent bound 

(3.47) |(L^- - Xr'd^Jkn^,) < Cp-i-i/^+n/l/lLi(xO + Ifk^i.,)], 
for all 2 < p < oo, (3 = 0,1, and e > 0. 

Proof. Following \Z2\ IZ3|. define the curves 

(|,A)(p,|,f) ■.= {pi,pif-9,p^). 
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where e M'^ G M and (^,f) G S'^: + |fp = 1. As {p,£,,f) range in the compact set 
[0, 6] X S'^, A) traces out the portion of the surface contained in the set + |Ap < 6. 
Thus, using and L°° estimates obtained in previous sections with 7 = and applying 
the interpolation inequahty between and L°° spaces, we obtain the proposition in the 
case /3 = 0. 

Now, recaUing that W = ^VZ and all coordinate transformation matrices are uniformly 
bounded, the refined bounds of Z therefore imply improved bounds for W and thus U. 
Bounds for L^, 2 < p < 00, are obtained by interpolation inequality between and L°°. 
Hence, we have proved the bounds for /? = as claimed. 

For /3 = 1, we expect that dx^f plays a role as "p/" forcing. Recall that the eigenvalue 
equations (L^ — X)U = dx^f read 

LoU 

{B'^Ux,)x^ - {A^U)x, -i ^'^JU + i ^ B^^^.Ux, 

(3.48) j^i j^i 

+ iY^{B^^i^U)x, - B^%CkU -XU = dxj. 
fe^i j,k^i 

Now modifying the nice argument of Kreiss-Kreiss presented in \KK\ IGMWZT] , we write 
U = V + Ui, where V satisfies 

(3.49) {Lo-X)V = dxJ, XI GM. 

Noting that and B^^ depend on xi only, we thus obtain by one-dimensional results (see 
|MaZ3t IZ3] ) the following pointwise bounds on Green kernel of X — Lq, 

(3.50) \dy,Gl{xuyi)\<Ce'P\^^-y^\{p + e-'\y^\). 
Hence, employing Hausdorff-Young's inequality, we obtain 

(3.51) + IK-Jl^xi) < Cp-'/P[p\f\Liix,) + 
for all 1 < p < cxD. 

Now from Ui = U — V and equations of U and V, we observe that Ui satisfies 

(3.52) {L^-X)Ui = L{V,Vx,), 

where L{V,Vx,) = pOi\V\ + \VxA). 

Therefore applying the result which we just proved for /3 = to the equations (|3.52p . we 
obtain 

\Ui\LPix,) < Cp-'-'/P+^[\LiV,Vx,)\mx^) + \LiV,Vx,)\L^^x,) 
(3-53) <Cp-'-'/P+^p[\Vk. + \VxAL. 

<Cp-l/P+^[|/|^i(,^)+p-VlL-(xi)]- 
Bounds on V and Ui clearly give our claimed bounds on U by triangle inequality: 

\U\lp < \V\lp + \Ui\lp. 

We obtain the proposition for the case (3 = 1, and thus complete the proof. 
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□ 



3.6. Estimates on the solution operator. In this subsection, we complete the proof of 
Proposition 12. li As mentioned earlier, it suffices to prove the bounds for Si{t), where the 
low frequency solution operator Si{t) is defined as 

(3.54) := pip e-^'^HL, - X^dX^ 

Proof of bounds on Si{t). Let u{xi,^, A) denote the solution of {L^—X)u = f, where f{xi,S,) 
denotes Fourier transform of /, and 

1 



u{x,t) ■.= Si{t)f 



(2'KiY J\l\<r Jr« 



Using Parseval's identity, Fubini's theorem, the triangle inequality, and Proposition 13.5 
we may estimate 



\u\ 



1 



2d 



e^^u{xi,^,X)dX 



re 



d^dxi 



< 



{2n) 



2d 



r« 



2 _ 



<c[\f\mx) + \fy,-? 



re 



^3leXtp-3/2+e^^ 



d^. 



Specifically, parametrizing F^ by 

A(e, k)=ik- 9i{k^ + ieT), ^ e M, 

we estimate 



r« 



2 _ 



di < 



-ei{k^ + \i\^)tp-3/2+e^j^ 



d^ 



d^ 



noting that f^d-i & is finite, provided a < d — 1. 

Similarly, we estimate 



""I r2,oo (t) < 



(27r) 



2d 



JkeXt 



r« 



u{xi,i, >^)\L^{xx)dX 



di 



re 



di 



where, parametrizing F^ as above, we have 



dC< J^e-"' 
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Finally, we estimate 

1 ' ' ^SReAti', 



where, parametrizing as above, we have 



< ^^-(d-l)/2-e/2_ 

The xi— derivative bounds follow similarly by using the version of the — > LP estimates 
for /3i = 1. The x— derivative bounds are straightforward by the fact that d^f = (ii)^ f- □ 

4. Two-dimensional case or cases with (H4) 

In this section, we give an immediate proof of Theorem 1 1.7i Notice that the only assump- 
tion we make here that differs from those in |NZ2j is the relaxed Hypothesis (H3'), treating 
the case of totally nonglancing characteristic roots, which is only involved in low-frequency 
estimates. That is to say, we only need to establish the — > bounds in low-frequency 
regimes for this new case. 

Proposition 4.1 (Low-frequency bounds; |NZ2j . Proposition 3.3). Under the hypotheses 
of Theorem for A G F^ (see (j3.46p ) and p := |(^,A)|, 9i sufficiently small, there holds 
the resolvent bound 

(4.1) \{L^ - A)-iaf < C72p-'/^[/|/|Li(.,) 
for all2 < p < oo, (3 = 0,1, and 72 is the diagonalization error (see |Z3| . (5.4O)) defined as 

(4.2) ^2 := 1 + ^ [p-i|9mA - rjf{0\ + p] ^'''~\ 
with V^,Sj as in (H4). 

Proof. We only need to treat the new case: the totally nonglancing blocks Qf. But this is 
already treated in our previous subsection. Subsection 13.4.41 yielding 

I |2 I 2| |2 ^ I zT" |2 

^^^^ +P \uHt+\L^ ^ \FHt+\L^, 

\uHt-\L°- +p'^\uHtJl2 < \uHt-{0)\'^ + \FHt Jli, 

where the boundary term |M_H't_(0)p can be treated by applying the stability estimate 
(j3.6p . Thus, together with a use of the standard interpolation inequality, we have obtained 

(4.4) \uHt\LPix^) < C72/3"V|l1(^'i)' 

for all 2 < p < 00 and 72 defined as in (j4.2p . yielding (j4.ip for f3 = 0. For (3 = 1, we can 
follow the Kreiss-Kreiss trick as done in the proof of Proposition 13.51 completing the proof 
of Proposition 14.11 □ 
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Proof of Theorem \1.7\ Proposition 14.11 is Proposition 3.3 in |NZ2j with an extension to the 
totally nonglancing cases. Thus, we can now follow word by word the proof in |NZ2j . 
yielding the theorem. □ 

Appendix A. Genericity of (H4') 

Genericity of our additional structural assumption (H4') is clear. Indeed, violation of the 
condition would require d equations: d^.Xk{(,) = for all j = 1, • • • , d, whereas only d — 1 
parameters in ^ G M'^ \ {0} are varied as ^ may be constrained in the unit sphere S'^ by 
homogeneity of A(^) in ^. 

Finally, we give the following counterexample of Kevin Zumbrun in the two-dimensional 
case for which the hypothesis (H4') fails. 



Counterexample A.l. Let 

(A.l) 



1 

1 



A. 




1 



Then both Ai and A2 are clearly symmetric and do not commute. However, at ^1 = 0, the 
matrix ^1^1 + ^2^2 has an eigenvalue (A(^) = 0) such that VA = 0, violating (H4'). 

Counterexamples for higher-dimensional cases can be constructed similarly. 
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